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ABSTRACT 
There  has  been  a  lot  of  activity  in  recent  years  on 
the  study  of  Lanches ter-type  combat  models,  especially  from 
a  deterministic  standpoint.   We  consider  some  of  these  models 
in  a  stochastic  framework  and  indicate  how  the  appropriate 
deterministic  model  can  be  recast  stochastically.   Techniques 
for  obtaining  the  corresponding  solutions  to  the  resultant 
differential-difference  equations  are  discussed.   These 
techniques  are  similar  to  those  developed  for  use  in  other 
population  process  modelling  situations  such  as  epidemic 
theory  and  competition  models.   These  are  used  to  give  the 
actual  solution  for  the  stochastic  model  of  the  original 
Lanchester  model. 


KEYWORDS:   Lanchester,  deterministic,  stochastic,  state 
probabilities,  duration  time. 


STOCHASTIC  LANCHESTER-TYPE  COMBAT  MODELS  I 

by 

L.  Billard 
Florida  State  University 

1 .   INTRODUCTION 

Since  World  War  II  there  has  been  much  activity  on  the 
study  of  combat  models  based  on  the  original  work  of  Lanchester 

(1914) .   The  so-called  Lanchester  models  describe  the  situation 
in  which  two  forces  are  in  combat  with  each  other  with  each  side 
losing  men  or  equipment  (tanks)  by  attrition  in  accordance  with 
some  preassigned  attrition  law.   These  laws  vary  from  the  very 
simplest  formulation  used  by  Lanchester  to  quite  complex  cases. 
An  extensive  review  of  these  different  situation  is  presented 
in  Taylor  (1978)  . 

Basically  the  model  for  the  combat  process  can  be 
described  in  terms  of  the  sizes  of  the  two  combat  forces.   Most 
of  the  work  in  the  literature  so  far  has  confined  its  attention 
to  studying  the  combat  process  in  a  deterministic  framework. 
While  such  an  approach  can  prove  useful  in  providing  broad  guide- 
lines as  to  the  behavior  of  a  given  combat  situation,  it  is 
likely  that  a  more  accurate  account  can  be  obtained  when  the 
process  is  viewed  stochastically.   Therefore,  in  this  paper  it 
is  shown  how  analogous  stochastic  combat  models  can  be  developed, 
and  the  appropriate  solutions  are  derived.   Specifically  we 
show  that  a  stochastic  analogue  of  the  Lanchester-type  combat 
models  is  nothing  but  a  particular  bivariate  death  process. 
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Furthermore  if  we  wish  to  extend  the  simplest  combat  models 
to  allow  for  reinforcements  we  will  have  a  bivariate  birth 
and  death  process. 

In  the  following  section,  we  describe  the  general  bivariate 
pure  death  process.   This  is  applied  to  the  particular  case  of 
the  original  Lanchester  model  in  Section  3.   In  Section  4,  we 

show  how  several  other  models  can  be  formulated  as  a  bivariate 
death  process.   These  are  models  which  generally  have  been 
considered  in  the  literature  from  a  deterministic  approach. 
Some  have  been  fully  solved  and  others  only  partially  solved 
(that  is,  for  some  finite  subset  of  the  governing  parameters) . 
We  note  that  using  the  results  of  our  Section  2,  formal  solu- 
tions can  be  found  to  the  stochastic  analogues  in  all  cases. 
Then,  in  Section  5,  we  consider  some  models  whose  stochastic 
formulation  is  that  of  a  bounded  bivariate  birth  and  death 
process.   Finally  we  consider  combat  duration  time  in  Section  6. 
Numerical  results  and  comparisons  for  varying  parameter  values 
and  battle  force  sizes  will  be  presented  in  a  companion  paper. 


2.   BIVARIATE  DEATH  MODEL 

Let  B(t)   and   R(t)   be  the  size  of  the  Blue  and  Red 
forces  at  time  t,  respectively,  and  let  particular  state  values 
be   b   and   r,  respectively.   Typically,  we  may  be  interested 
in  the  number  of  men,  or  of  tanks,  or  of  both  men  and  tanks, 
etc.   We  interpret  "size"  as  the  number  of  units  (men, 
tanks,  etc.)  under  discussion.    Suppose  that  at  time   t  =  0, 
B(0)  =  BQ   and   R(0)  =  RQ . 

The  two  forces  are  in  combat  with  each  other  and  units 
are  lost  according  to  various  attrition  laws,  the  particular 
attrition  law  depending  on  the  actual  modelling  situation  at 
hand.   For  general  attrition  rates   a(b,r,t)   and  3(b,r,t)  for  the 
Blue  and  Red  forces,  respectively,  the  general  deterministic 
model  is  such  that  the  following  equations  are  satisfied: 


(1) 


db      ,,    .  v 

-^  =  -a (b,r, t)  , 


||  =  -S(b,r,t)  . 


When  we  view  the  process  stochastically,  that  is,  when 
B(t)   and   R(t)   are  taken  to  be  random  variables,  the  attrition 
terms   a(b,r,t)   and  3(b,r,t)   translate  into  infinitesimal 
transition  probabilities  according  to  the  equations: 


P(B(t+h)  =  b-1,  R(t+h)  =  r|B(t)  =  b,  R(t)  =  r} 
=  a  (b,r, t)h  +  o(h)  , 

(2)      P{B(t+h)  =  b,  R(t+h)  =  r-l|B(t)  =  b,  R(t)  =  r} 
=  6(b,r,t)h  +  0(h), 

P{two  or  more  changes  in  (t,t+h)}  =  o(h)  , 

and  hence 

P{B(t+h)  =  b,  R(t+h)  =  r|B(t)  =  b,  R(t)  =  r} 
=  1  -  (a(b,r,t)  +  @(b,r,t)}h  +  o(h)  , 


where   lim,  n  o(h)/h  =  0  . 
If  we  write 


p.   (t)  =  P{B(t)  =  bf  R(t)  =  r}  , 
d  /  r 


then  the  forward  differential-difference  equation  governing  this 
process  is  (there  is  a  corresponding  backward  equation  also) 


(3)      A.  pb   (t)  =  -{a(b,r,t)  +  3(b,r,t)}  Pb#r<t) 


+  a(b+l,r,t)  pb+1  r(t)  +  3(b,r+l,t)  Pj,  r+1(t) 


for   (b,r)  6  A  =  {  (bf  r)  :  0  <_  b  £  BQ ,  0  <_  r  £  RQ } ,   and  where 


p,    (t)  =  0   whenever   (b,r)  /  A.   The  initial  conditions 
are 

%,R0(0)  =  l'   pb,r(0)  "  °   for   (b'r)  *    (VV' 


As  written  here  the  transition  probability  generators 
a(b,r,t)   and   3(b,r,t)   are  any  (positive)  function  of  the 
state  of  the  process  (bfr)  and  any  function  of  time  t.   Fre- 
quently, the  models  of  interest  are  such  that  these  generators 
are  time  independent,  that  is,  a(b,r,t)  =  a(b,r)   and 
3(b,r,t)  =  ,8(b,r).   When  this  holds,  the  set  of  equations  (3) 
is  just  the  bivariate  pure  death  process  described  by  Billard 

and  Kryscio  (19  77) .   They  then  give  the  solution  for   p,   (t) 

d  ,  r 

for  this  equation  for  any  generalized   a(b,r)   and   3(b,r) 
provided  that 


(4)    a(b1,r1)  +  60^,1^)  f   a(b2,r2)  +  3(b2,r2) 


for   (b-^r.^  t*  (b2,r2) 


Thus, 


(5) 


B0   R0 


Pb  r(t)     I         I      c,(m,w|B0,R  )  c2(m,w|b,r)  exp{d(m,w) t} , 
'        m=b  w=r 


where 

d(b,r)  =  -a(b,r)  -  3(b,r) 

and  where   c, (m,w|Bn,R_)  and   c?(m,w|b,r)   are  given  in  Billard 

and  Kryscio  (1977,  eqns .  10-11). 


Most  processes  of  practical  interest  are  such  that  (4) 
holds  naturally.   In  obtaining  their  solution,  Billard  and  Kryscio 
(1977)  essentially  exploit  and  utilize  a  more  general  theorem 
given  in  Severo  (196  9a) .   In  the  event  that  the  restriction  (4) 
does  not  hold,  the  solution  can  still  be  found  by  using  the  Severo 
theorem  directly.   Likewise,  when  the  transition  generators  are 
functions  of  time  t,  direct  use  of  Severo 's  theorem  yields  the 
required  solution  to  (2) .   In  any  event,  exploitation  of  the 
underlying  structure  which  expresses  itself  as  a  partitioning 
scheme  as  used  in  Billard  and  Kryscio  (1977)  and  Billard  (1980), 
has  the  advantageous  effect  of  reducing  considerably  the  degree 
of  complexity  that  a  first  glance  at  Severo' s  result  suggests 
is  involved. 


3.   LANCHESTER'S  MODEL 

The  simplest  formulation  for  combat  models  is  that  of 
Lanchester  (1914)  where  it  is  supposed  that  losses  on  each  side 
are  proportional  to  the  size  of  the  opposing  force.   This  is  the 
case  of  "aimed"  fire.   That  is, 

a(b,r,t)  =  a(b,r)  =  y^ 
and 

S(b,r,t)  =  3(b,r)  =  y2b  • 

For   convenience,   we    rescale    so    that      a(b,r)   =  Ar      and      3(b,r)    =b 
where   now      A      is    the    relative    effectiveness    of    the    Red    force    to 
the   Blue    force.      Then,    the    appropriate    differential-difference 
equation   governing    this   process    is    simply,    from    (3), 


<6)  &  Pb,r(t)    =    -(Xr  +  b>       Pb,r(t)    +    b     Pb,r+l(t)  +  te   Pb+l,r(t) 


for   (b,r)  d   A   and  where  we  note  that  for  boundary  values 
(b,r)   suitable  adjustment  is  necessary.   In  matrix  -Form,  we  mav 
write  (6)  as 

(|rPb,r(t))  =  B(pb;r(t)) 

where   B   is  the  matrix  of  coefficients  (of  transition  generators) 
Thus  in  the  particular  case  that   B_  =  4 ,  R.  =  2   and   A  =  .8, 
B   is  given  in  Table  1 . 

The  solution  to  (6)  is  given  by  (5)  where  it  is  readily 
verified  that 


(7a)   cx (mfw|BQ/R0) 


B  -m 
(Aw)  U   /(B0-m)l   , 


W  =  RQ, 


Rq-w+1 


V-w      il=1 


X2  X   °    R-1   °(Rn-l)  2      1- 

l 


1R0-w+l~1R0-w 
w 


(B0  "  m  +  1  ~  \-^  '■ 


RQ-w 

x  n 

v=l 


(Bo  +  x  -  V 

lBn-m+l-i  +A (R.-w+l-v) } 
0      v    0 


i  -1 
v 

x     n     {BQ-m+l-j  +A (RQ-w+l-v) } 


-1 


,   w  <  RQ  , 


with   i^.  =  1   and   i_    ,,  =  B.-m+l;  and 

U  Rjn.-W+l  U 


(7b)  c2(m,w|b,r) 


(-l)m-b(Aw)m-b/(m-b)l   , 


r  =  w  , 


i    .-,           i0                    ,  1,-1,,      i»-i,        l    ,,-i 
w-r+1      2      l    ,,  -i     ,  m-b   1   0  ,   ,  »  2   1        w-r+1   w-r 
y        S       (    1)    w~r+*-      w~r  ^    w      (w-1)      r 


|i      n=l        i   =1 

w-1  1 


{A  (w-r)  +  m-  b>  './{A  (w-r)  +i         =T7T 

w-r 


w-r 

x     n 

v=l 


(m-i    +1) 
v 

(1-vA-i    ) 


i    -1 
v 


V        1     =1        , 

Jv      v-1 


(A(l-v)-jv) 


-1 


r    <   w    , 


with      i  , ,    =   m-b+1 . 

w-r+1 


In  matrix  notation,  the  set  of  state  probabilities  may 
be  written  as 

(pb,r(t))  =  £  -(t)  ' 

where  the  vector   e(t)   has  elements   {exp(b.t)}   with   b. 
being  the  ith  diagonal  element  of   B   and  the  elements  of  the 
matrix  C   being  determined  from  (7) .   In  our  example,  the 
corresponding   C  matrix  is  given  in  Table  2.   Thus,  for  example, 


P(B(t)  =  2,  R(t)  =  2}  =  1.28e  5,6t-2.56e  4*6t  +  1.283e  3,6t 


Further  details  on  how  to  calculate  these  quantities  algorithmically 
will  be  discussed  in  the  forthcoming  paper  dealing  with  numerical 
calculations  and  comparisons  generally. 

We  notice  that  the  underlying  structure  for  both   c,  (  * ) 
and   c-f*)   is  very  similar,  and  furthermore  within  each  formula 
the  structure  for  each   i    term  is  also  similar.   This  allows 
computer  calculation  of  these  quantities  to  run  without  too  much 
difficulty.   This  effect  is  especially  apparent  in  the  particular 
case  that   X  =  1,  that  is,  when  the  two  forces  are  equally 
effective.   In  this  case  we  find  that 


c1(m,w|B(),R0)  •  (BQ    +    RQ    -   m   -   w)  '. 


B0-m 

w  , 


w   =    RQ    , 


R   -w+1 


R0-w=l 


x2         i,  -in  io~i-\ 

11=1 


RQ-w 

X         nn        (B0    +    l    ~    iv)     ' 
v=l 


and 


c2  (m,w|b,  r)  •  (m+w-b-r)  '. 


w 


1RQ-w+l    1RQ-w 


w    <    RQ    ; 


,    . »  m-b     m-b 
(-1)  w        , 


"w-r+1 

I     •■ 

Lw-r=1        il=1 


2  r.         i,-in  1.-1, 

v       /    ,x  m+w-b-r      1      0,      , v    2      1 
)_       (-1)  w  (w-1) 


r  =  w    , 


w-r+1     w-1 


w-r 
x      n      (m  +  1   -   i    )    , 

-i  v 

r=l 


r    <   w 


Gye    and  Lewis    (19  76)    have    argued   that   it    is    quite    reasonable    that 

in    the   Battle   of  Trafalgar,    A   =   1.      Other   such    situations   can  be 

easily   visualized. 

Once   expressions    for      p,        (t)    have  been    determined,    other 

quantities   of   interest   can   be   obtained.      Suppose   we    are    interested 

in   the   event   that   in   time    t   the   Blue    force    loses      x      units   while 

the    Red   force    loses   none.      Substitution   in    (5)    and    (7)    yields 

the    result 
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(8)    p{B(t)  =  BQ    -    X,  R(t)  =  RQ} 

=  PB  _XfR    (t)  =  exp[-(B0-x+RQA)t]  (AR0)X(l-e"t)X/x'.  , 


X  —   U  /   •  •  •  i  f\       * 


Likewise,  if  we  are  concerned  with  the  probability  that  Red 
loses  no  units  (regardless  of  the  number  of  loses  on  the 
Blue  side)  Up  to  time  t,  we  have 

B0 
(9)  P(R(t)  =  R  }  =  I      p        (t) 

U     x=0   B0  X'R0 

which  can  be  determined  from  (8). 

These  results  (8)  and  (9)  could  be  likened  to  the  model 
described  by  Gaver  (1979)  in  which  a  force  of  initial  size   B 


0 


attacks  a  bastion  or  stronghold  of  size   Rn .   Gaver  assumed 


that  the  Red's  stronghold  is  sufficient  to  guarantee  no  loss  of 
Red  units.   He  then  considers  the  size  of  the  Blue  force  at  time  t 
using  always  the  deterministic  approach.   Gaver  looks  at  two 
cases.   One  case  assumes  that  the  attrition  rate  of  the  Blue 
force  is  a  result  of  unaimed  fire  by  the  Red  force.   Therefore, 
in  the  terminology  of  the  present  paper,  the  (stochastic)  trans- 
ition rate  is 

a(b,r,t)  =  pu(R0/BQ)b, 
where   pu   is  the  attrition  parameter.   In  the  other  case,  Red 
employs  aimed  fire  and  hence  the  combat  is  modelled  so  that 

a  (b,  r , t)  =  p  b 
a 
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for  attrition  parameter   p  .   In  either  case,  Gaver's  process 

a 

is  slightly  different  from  the  Lanchester  results  (8)  and  (9) 
since,  although  the  quantity  of  interest  includes  the  fact  of  no 
loses  to  the  Red  force,  it  is  nevertheless  possible  probabilisti- 
cally that  the  Reds  do  lose  some  units  by  attrition.   In  Gaver's 
model  it  is  assumed  the  Reds  cannot  lose  any  units  whatsoever. 

Returning  to  the  Lanchester  model,  we  can  easily  obtain 
the  expected  number  of  units  lost  by  time  t  by  the  Blue  force 
given  that  the  Red  force  has  no  loses.   Let  this  random 
variable  be  denoted  by   X(t) .   Then, 


B0 
E(X(t)}  =  I      x  P{B(t)  =  B   -  x|R(t)  =  RQ} 
x=0 

Bo  Bo 


"   £   X  PB  -x  R  (t)/  E   PB  -x,R  (t) 

x=0    B0  X/Ro    x=0  *0  x'Ko 


=  AR0(l-e  t)    S(B0-1)/S(B0)  , 


where 


S(p)  =  I    {XRn(l-e"t)et}y/Y'.  • 
y=0 

If   B    is  sufficiently  large  so  that   S(BQ)  ~  S(BQ-1),  we  have 


E{X(t)}    ~    XR0(l-e    fc)     . 


He 


nee,    the  eventual   expected  number   of    loses    on    the    Blue    side    is 


lim  E{X(t)>    =    ARQ    . 
t  •*■  °° 


This  result  is  what  we  would  expect  intuitively. 
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We  can  relate  the  stochastic  and  deterministic  models 
through  the  expectations  since  the  attrition  rates  are  linear 
functions  of  the  variables  involved.   Thus,  if  we  multiply  (6) 
throughout  by  b   and  sum  over  all   (b,r)   values,  we  have 


A.  E(B(t) }  =  -AE{R(t)} 


Likewise,  multiplying  by   r   and  summing  over   (b,r)   gives 


^  E{R(t)}  =-  E{B(t)}  . 


Comparison  with  the  deterministic  equations  shows  that  the 
expected  values  of  the  stochastic  variables  equals  the  solution 
of  the  deterministic  equations.   We  note  this  correspondence 
does  not  necessarily  hold  when  the  attrition  rates  become  non- 
linear in   b   and/or   r. 
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4 .   OTHER  COMBAT  MODELS 

In  this  section  we  briefly  describe  how  other  combat  models 
can  be  expressed  stochastically.   There  have  been  many  models 
proposed  in  the  literature.  An  extensive  and  exhaustive  review 
of  such  models  has  been  made  by  Taylor  (19  78)  to  which  we  refer 
the  reader  for  further  details,  elaboration,  and  justification. 
These  models  have  by  and  large  been  studied  from  a  deterministic 
viewpoint  only.   We  confine  outselves  here  to  establishing  certain 
stochastic  analogues  of  some  of  these  processes.   Specifically, 
we  obtain  the  appropriate  format  for  the  transition  probability 
generators   a(b,r,t)   and  3(b,r,t). 

Once  the  generators  have  been  established,  it  is  then  a 
relatively  simple  and  direct  procedure  to  solve  the  corresponding 
differential-difference  equations  as  discussed  earlier  in  Section  2 
Obviously  the  simpler  the  form  for   a(b,r,t)   and  3(b,r,t),  the 
simpler  the  resulting  solutions  will  be,  such  as  we  saw  in  the 
previous  section  for  the  original  Lanchester  model.   Quite  clearly 
the  reduction  achieved  depends  on  the  actual  structure  of  these 
generators.   When   a(b,r,t)   and  S(b,r,t)   are  linear  functions 
of   b   or   r,  such  reduction  will  be  substantial  and  nice.   How- 
ever, even  the  most  complicated  cases  can  still  be  solved  since 
use  of  the  Severo  theorem  is  completely  general  without  any  con- 
fining restrictions  as  far  as  our  models  are  concerned. 

Taylor  and  Brown  (19  76)  establish  what  they  call 
"variable-coefficient  Lanchester-type  equations  of  modern  warfare" 
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where  the  attrition  rate  for  each  side  is  time-dependent.   By 
analogy,  we  find  the  transition  probability  generators  to  be 
given  by 


(10) 


a(b,r,t)  =  a(t)r    and    3(b,r,t)  =  b(t)b  . 


If   a(t)  =  Y-ih(t)   and  b(t)  =  y~h(t),  rescaling  of  time  from 
t   to   h(t)   produces  generators 


a (b,r, t)  =  Yxr 


and 


(b,r,t)  =  Y2b  ' 


which  is  the  Lanchester  model  of  Section  3. 

Range-dependent  attrition  rates  were  introduced  by 
Bonder  (196  7) ..   Here,  we  have 

a(b,r,t)  =  -a(d)r   and   3(b,r,t)  =  -3(d)b  , 

where   a(d)   and   3(d)   are  the  attrition  rates  for  the  Blue  and 
Red  forces,  respectively,  dependent  on  the  range  or  distance 
between  the  two  forces.   In  some  cases  of  tactical  interest, 
a (d)   and   3(d)   may  be  represented  by  the  function 


a(d)  = 


p,  (1  -  d/dn) 


0   , 


0  £  d  £  d2  ' 


d2  <  d  , 


and 
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3(d)  = 


p2d  -  a/d1)  z,    o  <_  d  <_  a1   , 


0   ,  d,  <  d  , 


where   d,   and   d~   are  the  maximum  effective  ranges  of  the  Blue 
and  Red  forces,  respectively,  and  where   p. ,  p2,   y.   and   y~   are 
appropriate  nonzero  constants.   If  at   t  =  0 ,  d  <_  M  in   (d- ,d2)  , 
that  is, when  combat  begins  both  forces  are  within  weaponary  range 
of  each  other,  and  if  there  is  no  retreat,  then  this  process  re- 
duces to  the  Lanchester  model  of  Section  3  where  now 

yi 
y±   =    Pi(l  -  d/d3_i)  \    i  =  1,2  . 

However,  it  is  equally  (perhaps  more)  reasonable  to  assume 
that  the  range  parameter   d   is  itself  time  dependent.   Bonder 
(1967)  takes 


d  e  d(t)  =  dQ  -  vt 


where   d.   is  the  initial  range  at   t  =  0   and  v   is  the  constant 
attack  (closure)  speed.   Thus,  the  attrition  rates   a(d)   and 
3 (d)   are  now  time  dependent. 

The  range-dependent  attrition  rates  of  Bonder  are  a 
particular  case  of  a  more  general  power  attrition  rate  class  of 
models  discussed  by  Taylor  and  Brown  (1976)  and  Taylor  and  Comstock 

(1977) .   This  general  class  is  such  that  the  transition  prob- 

pl 
ability  generators  assume  the  form   a(b,r,t)  =  k  (t  +  C)     and 
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^2 
8(b,r,t)  =  k2(t  +  C  +  A)   ,  where   A  >  0   is  called  the  offset 

parameter  and  C  ^  0  is  called  the  starting  parameter.  This 
accommodates  the  situation  in  which  the  forces  have  different 
maximum  effective  ranges  and/or  the  battle  begins  within  those 
maximum  ranges.  Taylor  and  Corns tock  (19  77)  give  the  solution 
to  the  deterministic  model  when  there  is  no  offset  that  is, 
A  =  0,  as  well  as  for  A  >  0  when  y,  =  y„  =  1  and  y,  =1, 
U~  =  2.  We  note  that  the  stochastic  solution  can  be  obtained 
for  all   u,   and   y?. 

A  distinguishing  characteristic  of  all  the  models  considered 
so  for  is  that  the  attrition  rate  of  each  force  is  proportional  to 
the  size  of  the  other  force  and  in  no  way  depends  on  its  own  size. 
It  is  quite  reasonable  to  expect  that  attrition  could  be  propor- 
tional to  the  size  of  both  opposing  forces.   This  is  especially 
so  far  unaimed  fire  (but  can  be  equally  argued  for  aimed  fire) 
when  it  is  realized  that  the  probability  that  a  unit  on  the  Red 
force,  say,  will  kill  one  unit  (out  of  the  remaining   B(t) 
units)  on  the  Blue  force,  increases  (decreases)  as  the  size   B(t) 
increases  (decreases)  since  there  are  more  (fewer)  units  on  which 
his  fire  may  fall.   In  this  context,  the  appropriate  transition 
probability  generators  in  a  time  independent  situation  would  be 
a(b,r,t)  =  -y-|br   and   3(b,r,t)  =  y  Jo>r .      We  note  that  here  the 
generators  are  time  independent  so  that  the  solution  (5)  holds 
directly.   However,  there  is  effectively  a  time  dependency  in 
that  the  actual  size  of  each  force  is  of  course  a  function  of  time. 


17 


This  model  is  mathematically  equivalent  to  the  Weiss  (1963) 
predator  prey  model  whose  solution  is  given  in  Billard  and 
Kryscio  (1977) . 

Let  us  now  consider  the  situation  in  which  one  or  both 
sides  has  a  supporting  system  which  itself  is  not  subject  to 
attrition.   Taylor  and  Parry  (19  75)  considered  a  particular  case 
of  such  a  process  which  when  presented  as  a  stochastic  process 
will  have  transition  probability  generators 

a(b,r,t)  =  -a(t)r  -  b1(t)b  , 
and 

3(b,r,t)  =  -b(t)b  -  a1(t)r  , 

where   a(t)   and  b(t)   have  the  same  interpretation  as  in  (10) 
and  where   a.  (t)  ^  0   and   b,  (t)  _>  0   represent  the  attrition 
coefficients  on  the  Red,  and  Blue,  forces  due  to  the  supporting 
fire  of  the  Blue,  and  Red,  forces,  respectively.   Other  situations 
come  readily  to  mind  such  as  a  model  in  which  the  attrition  from 
the  supporting  fire  is  also  dependent  on  the  size  of  its  target 
(opposing  force) .   Then,  we  would  have 


a(b,r,t)  =  -a(t)r  -  b^tjbr  , 


with   3(b,r,t)   similarly  defined  if  the  same  conditions  prevail 
for  the  Red  side. 
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Finally,  Helmbold  (1965)  concerned  himself  with  the 
situation  in  which  grossly  unequal  force  sizes  are  in  combat.   His 
modification  of  the  Lanchester  equations,  when  viewed  stochastically, 
gives  us 

a(b,r,t)  =  -a(t)  h(b/r)r 
and 

B(b,r,t)  =  -b(t)  h(r/b)b  , 

where   h(z)   is  interpreted  as  the  effectiveness-modification 
factor.   Helmbold  imposed  some  conditions  on  the   h(z).   However, 
we  note  that  for  the  stochastic  model  we  are  able  to  find  the 
appropriate  conditions  without  his  restrictions 

Quite  clearly  the  concept  of  supporting  fire  can  be  combined 
with  Helmbold1 s  model.   This  has  been  done  by  Taylor  (19  76)  in 
the  particular  case  that  h(z)  =  z    for  some  constant   c  >  0. 
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5.   BIRTH  AND  DEATH  MODELS 

One  feature  of  the  models  discussed  to  date  is  that  only 
losses  due  to  attrition  occur  to  the  combat  forces,  that  is,  they 
are  pure  death  processes.   However,  there  are  some  situations 
which  may  be  adequately  modelled  as  a  birth  and  death  process  as, 
for  example,  when  reinforcements  are  permitted.   If  the  maximum 
number  of  possible  reinforcements  is  known  in  advance,  then  we 
have  a  bounded  birth  and  death  process. 

As  before,  let  B(t)   and   R(t)   represent  the  size  of 
the  Blue  and  Red  forces  at  time  t,  respectively;  and  let   Bn  =  B(0) 
and   Rfi  =  R(0) .   Suppose  that  at  time  t  there  remain   B*(t)   and 
R*(t)   Blue  and  Red  reinforcements,  respectively  and  suppose 
B*(0)  =  B*   and   R*(0)  =  R* ,  that  is,   B*   and   R*   represent  the 
maximum  possible  reinforcements  available.   Thus,  after  time  t, 
B*  -  B*(t)   actual  reinforcements  have  been  added  to  the  Blue 
combat  force,  with   R*  -  R£(t)   similarly  defined  for  the  Red 
side.   Let   B*  +  BQ  =  B   and   R*  +  R   =  r.    Let 

X(t)  =  (B*(t),  B(t),  R*(t),  R(t))   with  realization   x =  (b*,b, r* , r) 
Then,  for  a  completely  general  model,  we  have  the  following 
infinitesimal  transition  probabilities  (corresponding  to  (2)) 


P{X(t+h)  =  (b*-l,b+l,r*,r) |x(t)  =  x)  =  A1(x,t)h  +  o(h)  , 

p{X(t+h)  =  (b*,b,r*-l,r+l) |x(t)  =  x}  =  A2(x,t)h  +  o(h)  , 

(11)   P{X(t+h)  =  (b*,b-l,r*,r) |x(t)  =  x}  =  U1(x,t)h  +  o(h)  , 

P{X(t+h)  =  (b*,b,r*,r-l) |x(t)  =  x}  =  y2(x,t)h  +  o(h)  , 


P{two  or  more  changes  in   (t,t+h)}  =  o(h)  , 
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and  hence 


P{X(t+h)    =   x|x(t)    =    x}    =   1    -      I      (xi(x,t)    +   yi(x/t)}h   +   o(h) 


i=l 


If  we  write 

p(x,t)    =    P{X(t)    =   x}, 

the   differential-difference   equation   governing    the    process    is 
(compare    ( 3) ) 


(12)         ^  p(x,t) 


=    "      I 


i=l 


{X±(x,t)    +  yi(x,t)  }    p(x,t) 


+  X1(b*+l,b-lfr*fr,t)  p (b*+l ,b-l ,r*, r , t) 
+  X2(b*,bfr*+l,r-lft)  p (b*,b, r*+l , r-1 , t) 
+   yi(b*,b+l,r*,r)    p(b*,b+l,r*,r,t) 

+   y2(b*,b,r*,r+l)    p (b* ,b, r* ,r+l , t)     , 


for 


x    €    8    = 


x:0    <_  b*    <_  B*,    0    <_   r*    <    R* 
0<b       <B,       0<r       <R 


and  where   p(x,  t)  =  0   whenever   x  /  B.   The  initial  conditions 
are 

p(B*,B0,R*/RQ,0)  =  1  ,   p(x,0)  =0    for   x  ^  (B0 'BQ ' R0 ' V  * 
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The  equation  (12)  can  be  solved  by  combining  the  techniques 
of  Severo  (1969a, 1969b) .  When  the  transition  generators  are  time 
independent,  that  is,   X.  (x,  t)  =  A.  (x)   and   y  .  (x,  t)  =  jj  .  (x)  , 
i  =  1,2,  an  explicit  solution  for  any  arbitrary   A. (x)   and 
y . (x) ,  i  =  1,2,  is  given  in  Billard  (1980). 

The   y.(x,t),  i  =  1,2,  terms  give  the  attrition  rates 
and  typically  we  would  expect  them  to  assume  the  same  forms  as 
the   a(b,r,t)   and   3(b,r,t),  respectively,  presented  earlier. 
The  terms   A. (x,  t),  i  =  1,2,  represent  the  rates  that  reinforce- 
ments are  added.   Though  written  generally  here,  it  is  reasonable 
to  expect  that  for  models  of  practical  interest,  these  generators 
will  be  a  function  of  the  size  of  the  combat  force  itself.   Thus, 
for  example,  the  case  where  more  reinforcements  are  brought  in 
when  the  combat  force  itself  becomes  small  could  be  reflected  by 
a  transition  generator 

Ax(x,t)  =  A(t)/b  , 
or  simply 

Ax(x,t)  =  A/b  . 

Recently,  Gaver  (1979)  investigated  the  problem  of  infor- 
mation flow.   In  his  model,  individual  units  start  off  with  un- 
aimed  fire  but  then  as  information  about  the  opposing  force  is 
gathered  the  fire  becomes  aimed  fire.   Thus,  in  the   X(t)   notation 
used  above,  B*(t)   is  the  number  of  Blue  units  still  using  unaimed 
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fire  at  time  t,  B(t)   is  the  number  of  Blue  units  using  aimed 
fire  at  time  t,  and   R*(t)   and   R(t)   are  the  corresponding 
quantities  for  the  Red  force.   If  all  units  improve  their  fire 
to  aimed  fire  before  they  themselves  are  killed,  we  have  exactly 
the  situation  modelled  in  (11)  and  (12) .   However,  it  is  more 
realistic  to  assume  that  some  of  the  units  never  learn  from  the 
information  gathered,  thus  remaining  with  unaimed  fire,  and  will 
themselves  be  lost  by  attrition.   To  model  this  situation,  the 
infinitesimal  transition  probabilities  (11)  still  apply  in 
addition  to 


and 


p{X(t+h)  =  (b*-l,b,r*,r) |X(t)  =  x}  =  y1(X't)n  +  °(h)  / 


P(x(t+h)  =  (b*,b,r*-l,r) |x(t)  =  x}  =  y2(x,t)h  +  o(h)  , 


and  where  now 


2 
P{X(t+h)  =  x|X(t)  =  x}  =  1  -   V  (a-  (x,t)  +  y  (x,t)  +  Y.  (x,t)  }h  +  o(h). 

i41  i  -     i-    'i- 


Thus,  for  each   i  =  1,2,  the  \.   (x, t)   represent  the 
information  flow,  the   y.(x, t)   represent  the  loss  of  aimed  units 
from  attrition,  and  the  y.   (x, t)   represent  the  loss  of  unaimed 
units  from  attrition.   Suitable  adjustment  of  equation  (12)  will 
give  us  the  appropriate  differential-difference  equation  for  the 
process.   This  latter  equation  an  then  be  solved  using  the  method 
of  Severo  (1969a, 1969b)  . 
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It  is  now  possible  to  establish  Gaver's  (1979)  model  for 
information  flow  under  mutual  attrition  in  a  stochastic  frame- 
work .   Thus ,  we  have 

X1(x,t)  =  auab*     and    A2(x,t)  =  3uar*  , 

where      a  and      3  represent   the    rate   of   conversion    from  an 

ua  ua  c 

unaimed  capacity   to   aimed  capacity   for   the   Blue    and   Red   forces, 
respectively; 

y1(x/t)    =    puar*b/B   +    paarb/(b   +   b*) 
and 

y2(x,t)    =   Yuab*r/R   +   yaabr/(r   +    r*)  , 

where      p  represents    the   attrition   rate   of   the   unaimed  Red 

^ua  ^ 

forces   against  the   aimed  Blue    forces,      p  represents    the   attrition 

aa 

rate  of  the  aimed  Red  forces  against  the  aimed  Blue  forces,  and 

Y     and  y  are  the  corresponding  attrition  rates  of  the  Blue 

ua        aa 

forces  against  the  Red  forces ;  and 

Y,(x,t)  =  p   r*b*/B  +  p   rb*/(b  +  b*) 
' 1  — '      Huu    /     Mau   ' 

v0(x,t)  =  Y   b*r*/R  +  y   br*/(r  +  r*)  , 
'  2  —  ' uu  'au 


and 


where   p     and   p     represent  the  rate  of  attrition  to  the  unaimed 
^uu       Kau    ^ 

Blue  forces  due  to  the  unaimed  and  aimed  fire,  respectively,  from 
the  Red  forces,  and   y    and  y  are  similarly  defined  attrition 

rates  of  the  Red  forces. 

24 


6 .   COMBAT  DURATION  TIME 

The  previous  sections  have  been  concerned  primarily  with 
techniques  for  obtaining  the  state  probabilities  of  the  underlying 
distribution.   A  knowledge  of  these  probabilities  allows  most 
other  quantities  of  interest  to  be  derived.   One  such  quantity 
is  the  combat  duration  time  or  the  expected  length  of  time  for 
the  battle.   While  this  can  certainly  be  obtained  from  our  earlier 
results,  we  present  here  a  relatively  easy  straightforward  but 
recursive  method  for  its  derivation. 

Suppose  that  at   t  =  0,  the  Blue  force  has   B_   units 
and  the  red  force  has   Rn   units.   Let  us  assume  that  the  combat 
ends  when  the  Blue  force  is  reduced  to  a  size  of   B'   units  and/or 
the  Red  force  is  reduced  to  a  size  of   R'   units.   Hence,  the 
permitted  state  space  for   (b,r)   is  now  A'  =  {(b,r)  :B'  <_   b    B^ , 
R1  <_  r   <_   Rfl}   instead  of  the   A   used  in  the  previous  section. 
(The  adjustment  from  A   to   A '   in  the  results  given  earlier 
is  trivial.)   Let  T(b,r)   be  the  expected  combat  duration  time 
from  the  time  that  there  are   b   Blue  units  and   r   Red  units 
remaining.   Then,  T(Bn,R_)   will  be  the  overall  combat  duration 
time.   For  ease  of  illustration,  let  us  suppose  we  have  the  simple 
Lanchester  combat  model  where   a(b,r,t)  =  Ar   and   3(b,r,t)  =  b. 
More  general  models  are  treated  analogously. 

We  first  recall  that  the  basic  underlying  process  is 
Markovian.   From  (2),  and  basic  properties  of  Markov  processes, 
we  see  that  when  the  system  moves  from  the  state   (b,r)   it  moves 
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to  the  state   (b-l,r)   with  probability   Xr/ (b  +  Xr)   or  to 
the  state   (b,r-l)   with  probability  b/(b  +  Xr) .   Finally,  before 
leaving  the  state   (b,r)   it  will  have  stayed  there  for  an  average 
time  of   l/(b  +  Xr) .   Therefore,  we  can  write 

(13)    T(b,r)  =  l/(b  +  Xr)  +  Xr/(b  +  Xr)  T(b-l,r)  +  b/ (b +  Xr)  T(b,r-1), 

for  (b,r)  €  A'.  In  (13),  T(b,r)  =0  for  (b,r)  /  A'.  Clearly 
then,  by  starting  at  (b,r)  =  (B',R' )  and  proceeding  with 
(B',r),  r  =  R'+l,  ...,  RQ ,  returning  to  (B1 +1, r) ,  r  =  R* , . . . ,RQ , 
and  so  on  to  (BQ/r) ,  r  =  R ',..., R^,  we  can  find  T(Bn,Rn).  Note 
that  in  fact  we  have  a  matrix  of  T(b,r)  values  which  give  us 
the  duration  time  from  any  intermediate  stage  (b,r)  in  addition 
to  the  overall  duration   T(B0,Rn). 
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